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ABSTRACT 

We use the lSternberg et alJ (120141) theory for interstellar atomic to molecular hydrogen (Hl-to-H2) conversion 
to analyze Hl-to-H2 transitions in five (low-mass) star-forming and dark regions in the Perseus molecular 
cloud, Bl, BIE, B5, IC348, and NGC1333. The observed HI mass surface densities of 6.3 to 9.2 M© pc“^ 
are consistent with Hl-to-H2 transitions dominated by Hl-dust shielding in predominantly atomic envelopes. 
Eor each source, we constrain the dimensionless parameter aG, and the ratio luy/n, of the EUV intensity to 
hydrogen gas density. We find aG values from 5.0 to 26.1, implying characteristic atomic hydrogen densities 
11.8 to 1.8 cm“^, for Ijjy ^ 1 appropriate for Perseus. Our analysis implies that the dusty HI shielding layers are 
probably multiphased, with thermally unstable UNM gas in addition to cold CNM within the 21 cm kinematic 
radius. 

Subject headings: ISM individual objects (Perseus) - ISMxlouds - photon dominated regions (PDR) - galax¬ 
ies: star formation 


1. INTRODUCTION 

Conversion of hydrogen gas from atomic (HI) to molecular 
(H2) form is of critical importance for the evolution of the 
interstellar medium (I SM) and for star-formation in galaxies 

(ISternberg et al.ll2Q14 , hereaf ter S14). _ 

Recently, iLee et al.l (1201 2h and iLee et al.l (1201 51 hereafter 
LI 2/Ll 5) analyzed Hl-to-H2 transitions in several subregions 
within the well-studied Perseus molecular cloud. The Perseus 
cloud is itself part of the nearby Taurus-Auriga-Perseus com¬ 
plex located at a distance of ^300 pc and is em bedded within 
the Per OB2 HI supershell (iBallv et al.ll2008h . The Perseus 
cloud mass is a few lO'^Mr:), and this includes th e HI plus H2 
(iBachiller & Cernicharolll98^ iKirk et al]l2006h . The cloud 
consists of an extended dusty HI envelope surrounding several 
condensations of dense Ho gas (ISancisi et al.iri974 ISargentI 

[T?7^ it Jngerechts & ThaddeiisIM 9871: iRidge et al.ll2nn6h . The 

overall angular size depends on the tracers used, and a char- 
acteristic diam^^eter b ased on 21 cm kinematics is ^ 80 pc 
(llmara & Blitzll20ll1) . 

LI 2/Ll 5 used the HI data provided by the Galactic 
Areci bo L-band Leed Array HI Survey ((TALLA-HI sur¬ 
vey; iPeek et al"]l2011h . together with far-infrared data from 
the Improved Reprocessing of th e IRAS Survey (IRIS; 
iMiville-Deschenes & LagacMl2005h and the V- band extinc¬ 
tion i mage provided by the COMPLETE Survey (iRidge et al.l 
120061) . to derive HI and H2 surface densities (Zhi and ZH2) 
on ^ 0.4 pc scales, for several hundred sight-lines towards 
five dark and (low-mass) star-forming regions within Perseus. 
These are Bl, BIE, B5 (dark), IC348, and NGC1333 (star¬ 
forming). Here we consider the data presented by LI5 for 
which the inferred HI column densities are corrected for 
(small, up to 20%) 21 cm optical depth effects. 

LI 2/Ll 5 a nalyzed the Per s eus d ata using the formalism 
presented by iKrumholz et al.l (I2009L hereafter KMT) for in¬ 
terstellar Hl-to-H2 transitions in optically thick media. In 
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this paper we use the simpler and more general theory pre¬ 
sented by S14 to reanalyze the Perseus observations. In § 1, 
we briefiy summarize the relevant S14 formalism. In § 2 we 
present and fit the LI5 observations of the Zh 2 /^hi ratios in 
Perseus. In § 3 we use the observed maximal HI mass sur¬ 
face densities towards each H2 cloud to constrain the control¬ 
ling dimensionless parameter aG and the characteristic HI gas 
densities in the atomic envelopes. The observed HI mass sur¬ 
face densities are consistent with Hl-to-H2 transitions domi¬ 
nated by Hl-dust shielding. The relatively low gas densities 
we infer suggests that the HI shielding layers are probably 
multiphased and are not pure CNM. 


2. THEORY 

2.1. HI Column Density 

S14 presented a general analytic formula for the steady state 
column density of photodissociated HI gas in optically thick 
clouds illuminated by LUV radiation, derived for planar ge¬ 
ometry and uniform density gas. Lor irradiation by isotropic 
fields, the total HI column density is given by 


Am = 2 X 


(m) 


In 


1 aG 


■ + i 


( 1 ) 


In this expression, the factor of two is for two-sided illumina¬ 
tion, (/i) = 0.8 is a geometrical factor for isotropic radiation, 
Gg is the dust-grain absorption cross section per hydrogen nu¬ 
cleon for 912-1108 A Lyman-Werner (LW) band radiation, 
and aG is the basic dimensionless parameter. As in S14 we 
assume that 


= 1.9 X 10“^^ (pgZg cm^ , (2) 


where is the grain abundance relative to standard Galactic 
ISM grain abundances. Thus, is the “dust-metallicity”. In 
Equation ^ (j)g is 3. factor of order unity depending on the 
specific grain composition, and the scattering and absorption 
properties. The dimensionless parameter 


aG = 


DoG 
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Jdiss 
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Fig. 1.— = Eh 2 /I^hi as a function of Ztot = ^hi +^H 2 for five regions: Bl, BIE, B5, IC348 and NGC1333. The red curves are the “best-fits” to the 

data. The median errors (not displayed) are ±15% and ±35% in Etot and respectively. The best-fitting values of Ehi are also indicated. 


TABLE 1 

Total HI surface densities. 


Cloud 

Ehi (M© pc“2 ) 

Nm (lO^^ cm-2) 

Bl 

lA 

9.3 

BIE 

12 

9.0 

B5 

6.3 

7.9 

IC348 

7.1 

8.9 

NGC1333 

9.2 

11.6 


where Dq is the free-space H 2 dissociation rate (s“^), G is 
the average H 2 self-shielding factor, R is the H 2 formation 
rate coefficient (cm^ s“^), n is the total hydrogen gas den¬ 
sity (cm“^), /diss = 0.12 is the mean dissociation probability 
per H 2 -absorbed LW band photon, Fq is the free-space LW 
photon fiux (cm“^ s“^), and w = 1/[1 + (2.640gZ^)^/^] is the 
normalized H 2 -dust-limited dissociation bandwidth (see S14 
for a detailed discussion of all these quantities). Physically, 
aG is the ratio of the Hl-dust absorption rate of the effective 
unattenuated H 2 dissociation fiux, to the H 2 formation rate. 

For H 2 formation on dust grains R = 3 x 10“cm^ s“^ 
and the dimensionless parameter aG can be expressed as 

= 1 C4 Uv _ 

- n/lOOcm-3 l + (2.64(/.^Zpi/2 ’ 

where / py is the field intensity relative to the mean 
iDrainel (llQVSh interstellar field, such that Fq = 2.07 x 
10^/uv cm“^ s“^ and Dq = 5.Sx 10“^^/uv s“^. The value 
of aG determines the nature of the Hl-to-H 2 transition and 
the size of the integrated HI column. In the “weak-field” 
limit aG C 1, A^hi is small, and the Hl-to-H 2 transition is 
controlled by H 2 -line and H 2 -dust absorption. In the “strong- 
field” limit aG ^ 1, A^hi is large, and the Hl-to-H 2 transition 
is dominated by Hl-dust absorption. Importantly, for optically 
thick clouds the total HI column density, A^he depends only on 
aG and Og. 


Equation ([T]) for the HI column density is for steady-state 
conditions such that the local H 2 destruction rate equals the 
formation rate, at every location. The equilibrium time-scale 
for HI /H 2 formation-destruction is 


1 

“ D + 2Rn ’ 


(5) 


where D is the local (attenuated) photodissociation rate. 
For molecular gas D/{2Rn) <C 1 and 4q \/{2Rn) ^ 5 x 
10^/n yr. For atomic gas D/{2Rn) 1 and 4q 1/Z). In 

free-space D = Dq and 4q ^ 5.5 x 10^//uv yf- 


2.3. Multiphased Gas 

For a multiphased CNMAVNM mixture of HI gas in which 
the heating is dominated by photoelectric emission from dust 
grains, the field intensity /uy and the densit y ncNM of the 
CNM are correlated, with (IWolfire et al.ll2003h 


'^CNM = 22. 



4.1 

TTsrzp^ 



cm 


( 6 ) 

Here Z^ q is the gas phase carbon-oxygen abundance rela¬ 


tive to the abundances at solar metallicity. That is, Z[. q is the 
“gas-phase metallicity”. In Equation (|6l), 0 cnm expresses the 
range for which the CNM can be in pressure equilibrium with 
the WNM. Typically, 0 cnm ^ 3. For multiphase conditions 
the WNM density is nwNM ^ O.OlncNM, and the gas is ther¬ 
mally unstable ( UNM) for densities between ^wnm and ^cnm 
( see Figure 9 of IWolfire et al.ll2003h . 

It follows from Equations (01) and ([6l) that for Hl-to-H 2 tran¬ 
sitions occurring in pure CNM at thermal pressures allowing 
multiphase conditions (KMT, SI4) 


, 1 ^ ^ 1 7'0.365 . . 7 ' . 

aG=(aG)cNM=2.58(--) (^) 


X 


/ 3 w 2.62 

^0CNM ^ 1 + (2.64(^gZ^) 



2.2. Time-Scale 
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In § 4 we will consider perturbations to {ocG)cnm by varying 
independently the parameters Z^, and Z^q, all for (^cnm 
in the realistic range of 2 to 5. 


3. I:h2/5:hi IN PERSEUS 

Equation ([T]) can be reexpressed as an HI mass surface den¬ 
sity 


Zhi = 6.71 



Mopc , 


( 8 ) 


where cr ^-21 = Og cm^). (In Equation ^ the con¬ 
tribution of helium to the mass is not included). If Ztot = 
Zhi + Zh 2 is the total hydrogen gas surface density, and ^h 2 = 
Zh 2 /^hi is the molecular-to-atomic mass ratio, then 


'^H2(^tot) = - 1 • (9) 


Eor optically thick clouds, ZHi(crg, aG) (as given by Equation 
El) is independent of Zfj 2 9 and t^*H 2 varies linearly with Z^q^ 
with slope 1 ^Hi- 

In Eigure[T]we plot the LI5 data for <^h 2 versus Ztot for 
the five dark and star-forming regions in Perseus Bl, BIE, 
B5, IC348 and NGC1333. A complete discussion of the data 
extraction methodology is presented in L12/L15. Eor each 
region, each data point corresponds to a distinct sight-line 
through the complex. Many of the sight-lines probe substan¬ 
tial columns of H 2 and pass well through the Hl-to-H 2 tran¬ 
sition layers. Eor example, ^h 2 approaches 10 in IC348 and 
NGC1333. 

We fit our Equation (|9]) to the data points using a standard 
weighted-least-squares procedure, and find the best-fitting to¬ 
tal HI surface density for each region. The total HI sur¬ 
face mass densities lie within the narrow range of 6.3 to 9.2 
M© pc“^, and are listed in Table [T] 

The (red) curves in Eigure \T\ are our best fits for <^h 2 ver¬ 
sus Ztot as given by Equation (|9]) and the HI surface densi¬ 
ties in Table [B The theoretical curves are in excellent agree¬ 
ment with the data. This implies that the sight-lines are in¬ 
deed probing optically thick complexes with complete Hl-to- 
H 2 transitions. 

Eor characteristic CNM, aG = {aG)cNM ~ 2.58 as given 
by Equation (|7]). Eor a standard cr ^-21 = 1.9 the HI column 
for pure CNM shielding is then Z^^^ = 4.0 M© pc“^, sig¬ 
nificantly smaller than the observed total HI columns. This 
implies that if the entire HI columns are contributing to the 
shielding of the H 2 cores in Perseus, these shielding columns 
must be multiphased, not just CNM. Alternatively, if the 
shielding is assumed to be entirely CNM, not all of the ob¬ 
served HI contributes to the shielding. We discuss these con¬ 
clusions in more detail in ^4. 


4. ANALYSIS 

4.1. Shielding and HI Gas Densities 

According to Equation ffl or ([8]) the total HI surface den¬ 
sity depends on just Gg and aG, so that curves of constant HI 
surface density may be drawn in the Gg versus aG parameter 
space. In Eigure [2] (all three panels) we plot the locus curves 
(in color) for the Zhi inferred for each of the five Perseus re¬ 
gions. As expected, for any Zhi a large aG requires a large 
Gg, and vice versa. 

The two horizontal dashed lines in Eigure [2] represent the 
range of grain absorption cross sections we consider in our 


analysis, from the standard Gg-2\ = 1.9 to a larger Gg-2\ = 
3.8. We consider an enhanced Gg because LI5 found that 
the visual extinction per hydrogen gas column in Perseus is 
Ay/Mil = 1.0 X 10“^^ mag cm^, about a factor of 2 larger 
than for standard Galactic extinction, and this may imply a 
correspondingly larger than usual dust-grain absorption cross 
section. A larger Gg could be due to (a) altered dust properties 
at a normal dust-to-gas mass ratio, i.e. = 2 and Z^ = 1, or 
(b) typical (diffuse) ISM dust but with a higher abundance, 
i.e (l>g = I and Z^ = 2. In any case, the gas-phase metallic- 
ity, Z'^ Q, appears close to Solar in Perseus (iHernandez et al.l 

Eor any assumed Gg the implied aG for each source may 
be read off the plots in Eigure [2l In Table [2| we list the range 
of inferred aG parameters for each region, for Gg- 2 \ = 1.9 to 
3.8. The aG are large (^1) and this implies that the attenu¬ 
ation of the photodissociating LW radiation is in the strong- 
field limit with Hl-to-H 2 transitions dominated by dust ab¬ 
sorption within the outer atomic envelopes (“HI -dust”). This 
as opposed to Hl-to-H 2 transitions controlled by H 2 -line self¬ 
shielding. 

Eor any given aG, the effective gas densities, n, in the HI 
gas depends on the assumed EUV radiation intensity Ijjy (see 
Eq.[4]). By “effective” we mean for a uniform density medium. 
As discussed by L12 within most of the Perseus system, the 
photodissociating radiation is dominated by the background 
Galactic light. This is consistent with the overall thermal in¬ 
frared dust emission temperatures (16 — 22 K) as well as with 
the anomalous microwave emissions (iTibbs et al.ll2011h . Eor 
EUV dust heating, Ijjy ^ 1 within a factor of 2. The radiation 
fields near IC345 and NGC1333 may be locally enhanced by 
the presence of one or two B5 V-type stars. 

In Table[2]we list the inferred gas densities assuming /uv = 
1, for Gg- 2 i = 1.9 and 3.8. The inferred densities scale lin¬ 
early with the assumed I\jy. Eor Gg -21 =3.8 the densities 
depend on whether (a) (l>g = 2 and Z^ = 1, or (b) = 1 and 

Zg = 2. The densities are a factor-two smaller for the sec¬ 
ond option (see again Eq. (H). Overall the effective HI den¬ 
sities range from ^ 2 to 10 cm“^. The gas densities in the 
H 2 cores are likely larger, enabling molecule formation on a 
time scales 1 / {2Rn) 5 x 10^/^ yr (see Equation ([5])), within 

the lifetime of the Perseus cloud ^ 10 — 100 Myr. In the HI 
layers the equilibrium time-scales are 1/D l /{2Rn) and a 
photodissociation steady state is achieved. 

In Table[2]we also list the range of length scales, i = Mn/^, 
for the atomic shielding envelopes given the observed HI 
columns. With exception of NGC 1333 (which may be in¬ 
fluenced by a B5V star) the derived length scales are com¬ 
parable, and within factors 2-3, with the overall ^ 80 pc HI 
kinematic size scale of the Perseus complex. The inferred 
sizes are perhaps too large for “option-b” (Z^ = 2, = 1) and 

more consistent with “option-a” (Z^ = 1, = 2) for which a 

larger dust cross section refiects an intrinsic variation in grain 
properties. 


4.2. Is the HI Multiphased? 

The gas densities that we have inferred above are lower than 
the densities expected for pure CNM as given by Equation ([b]), 
and are intermediate between ^cnm and ^wnm. This suggests 
that the observed HI columns are multiphased mixtures. If 
most of the HI extends to just the kinematic diameter of ^ 
80 pc (llmara & Blitzll201 1 I) much of the HI must be thermally 
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Fig. 2.— The observed HI contours in the Gg — aG parameter space. The horizontal dashed lines are for C7^-21 = 1.9 and 3.8. The grey strips are where 
aG = (ofG)cNM (Equation (7]) for 0 cnm in the range 2 — 5, (with 0 cnm increasing from right to left across the strip). The dust cross section C7^-2i = l-90gZ^ 
varies with the dust abundance and the intrinsic dust absorption properties (j)g. The left panel is for variations in (j)g assuming Z' = q = 1. The middle and 
right panels are for variations in Z^ with = 1, assuming Z^ q = Z' (middle) or Z^ q = 1 (right), see §4 for details. 


TABLE 2 

aG, VOLUME DENSITY RANGES, AND LENGTH SCALES, EOR C7g-21 IN THE RANGE 1.9-3.8 EOR/uv = 1- 


Source 

aG 

'i>s = i 

— 2 with Z' = 1 

Z' — 1 

— 2 with = 1 

n (cm“^) 

Length Scale (pc) 

n (cm“^) 

Length Scale (pc) 

B1 

6.5-26.1 

9.1-3.6 

33-84 

9.1-1.8 

33 - 169 

BIE 

6.1-23.8 

9.6-3.9 

30-74 

9.6-2.0 

30 - 149 

B5 

5.0-17.7 

11.8-5.3 

22-49 

11.8-2.6 

22-97 

IC345 

6.0-23.2 

9.8-4.0 

30-72 

9.8-2.0 

30 - 144 

NGC1333 

9.5-47.0 

6.2-2.0 

61-189 

6.2-1.0 

61 - 379 


unstable, possibly in a cooling transition from the WNM to 
CNM phases. 

These conclusions are also indicated by the positions of the 
grey strips in Figure [2l The strips show the regions in the pa¬ 
rameter space for which aG= {ocG) cnm for different (realis¬ 
tic) perturbations of the controlling quantities in Equation (|7]). 
In the left-hand panel the grey strip is for variations in (j)g 
from 0.5 to 3 assuming Z' = Zq q = 1. The middle panel 
is for metallicities between 0.1 and 3 assuming = ^co 
and with (j)g = 1. The grey strip in the right-hand panel is 
for variations in just Z^ from 0.5 to 3, but with Z^q = I 
and again = 1. The width of each strip corresponds to 
the range 2 < 0 cnm < 5, for CNM at multiphased conditions, 
with 0CNM increasing from right to left across the strips. 

It is evident from Figure [2| that for cjg -21 between 1.9 and 
3.8, the contours for the observed HI column densities are 
to the right of the grey strips, with aG always significantly 
greater than (aG)cNM (i-e., n < ^cnm) fof all types of per¬ 
turbations in the parameters Z^ and Z^q shown in the 
three panels. This implies that the HI cannot consist of pure 
CNM gas. For the HI to be fully CNM, the dust absorption 
cross section would have to be lower than expected, e.g. if the 
metallicity were reduced (see middle panel). 

Importantly, our conclusion that the HI cannot be pure 
CNM depends on the inclusion of the denominator w = 1 / [1 + 
{2M(j)gZ'gy^^] in Equation Q. As discussed by S14, w ac¬ 
counts for the reduction of the effective dissociation band¬ 
width by H 2 -dust. (This factor is not included in the KMT 


approximations.) For example, w = 0.4 for cjg_ 2 i = 1.9 with 
= 1 and Z' = 1. Excluding w in Equation dH) would shift 
the “CNM strips” to the right in Eigure [2] much closer to the 
HI contours. Without the H 2 -dust absorption term the inferred 
effective HI volume densities would be more than twice larger 
than listed in Table [2l and much closer to ^cnm as given by 
Equation (|6]) for /uv = 1 • 

5. DISCUSSION AND SUMMARY 

LI 2/Ll 5 used the “ spherical cloud” model de veloped by 
KMT (and updated bv lMcKee & Krumholzll201Qh to analyze 
the Hl-to-H 2 transitions in the various Perseus regions (e.g., 
see Eig. 11 in L15). In their analysis, each region consists 
of many individual spheres with H 2 cores surrounded by HI 
shells, and each sight-line probes the area-averaged mass ra¬ 
tios = (^H 2 )/ (^hi) as functions of (Etot) for each sphere. 
Eurthermore, LI 2/Ll 5 adopted the KMT ansatz that the HI 
shielding envelopes are dominated by CNM, and then esti¬ 
mated 0CNM for each region assuming n = ^cnm- 

As discussed in detail by S14 the differences in the pre¬ 
dicted HI columns and H 2 mass fractions for individual clouds 
are very small for plane-parallel versus spherical geometries. 
Ho wever, the Zkmt p arameter differs from the S14 aG (see 
also ISternberal 198 8li . since it does not include the H 2 -dust 
absorption factor w = 1 / [1 + (2.640gZ^) . Eor a given aG 

(Equation (H) and an assumed /uv the implied gas density n is 
increased if w is excluded. Thus, including the H 2 -dust term 
in aG is essential for determining how close the gas density 
is to the CNM density for multiphased conditions. 
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In their analysis, LI 2/Ll 5 assumed that = Z^q = 1. 
They also set (Jg-21 = LO corresponding to = 1 /1.9 in our 
Equation (|2]). It is evident from the left-hand panel of our Fig¬ 
ure [2] that for (Jg-21 = 1-0, the contours for the observed HI 
columns imply aG of 2.0 to 3.4, and are to the right of the 
grey CNM strip. Enforcing CNM for the HI would then re¬ 
quire (^CNM ^ 2. LI 2/Ll 5 derived larger 0 cnm (^5 — 10) for 
the various regions because the H 2 -dust term w is not included 
in the Zkmt they used. Without this factor the grey CNM 
strips in Figure [2| are shifted to the right. For (7^-21 = LO the 
inferred dcNU factors are then increased to the larger values 
found by L12/L15. 

In our analysis, we do not assume a priori that aG = 
(aG)cNM, instead we infer the effective HI gas densities. Our 
inferred densities range from ^ 2 to 10 cm“^ depending on 
the precise FUV intensity in Perseus, and on the assumed 
metallicities and FUV grain absorption cross section. These 
densities suggest that the HI shielding envelopes in Perseus 
are likely multiphased mixtures. If most of the HI is lim¬ 
ited to the kinematic 21 cm radius and shields the H 2 cores, 
then a significant fraction must be thermally unstable. Alter¬ 
natively some of the HI could be very extended WNM and 
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